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Abstract 

Supergiant Fast X-ray Transients (SFXTs) are HMXBs with OB supergiant companions. I review the results of the Swift 
SFXT Project, which since 2007 has been exploiting Swift's capabilities in a systematic study of SFXTs and supergiant 
X-ray binaries (SGXBs) by combining follow-ups of outbursts, when detailed broad-band spectroscopy is possible, with 
long-term monitoring campaigns, when the out-of-outburst fainter states can be observed. This strategy has led us to 
measure their duty cycles as a function of luminosity, to extract their differential luminosity distributions in the soft 
X-ray domain, and to compare, with unprecedented detail, the X-ray variability in these different classes of sources. 
I also discuss the “seventh year crisis”, the challenges that the recent Swift observations are making to the prevailing 
models attempting to explain the SFXT behaviour. 
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1. Introduction 

Supergiant X-ray binaries (SGXBs) hosting an accret¬ 
ing neutron star and an OB supergiant companion, are 
divided into classical systems, showing a strong X-ray vari¬ 
ability with an X-ray luminosity dynamic range of 10-50, 
and supergiant fast X-ray transients (SFXTs, Smith et al. 
2004; Sguera et al. 2005; Negueruela et al. 2006). The 
dozen or so members of the latter class (see Romano et al. 
2014c, for a recent review), have a quiescent luminosity 
of 10^^ erg s“^ (in’t Zand 2005; Bozzo et al. 2010) 
and display X-ray flares reaching 10^®-10^^ erg s“^. They 
are therefore identified based on their characteristic high 
dynamic range in X-ray luminosity, which reaches up to 
'^10^-10® times the range observed in classical systems 
(Sguera et al. 2005; Romano et al. 2015), even though 
the supergiant stars in SFXTs and classical SGXBs share 
similar orbital periods and spectroscopic properties. The 
origin of this different behaviour is still a matter of debate 
(see, e.g. Bozzo et al. 2013, 2015). Viable models involve 
extremely dense inhomogeneities (“clumps”) in the winds 
of the SFXT supergiant companions, compared to classical 
systems (in’t Zand 2005; Negueruela et al. 2008), the pres¬ 
ence of magnetic/centrifugal gates generated by the slower 
rotational velocities and higher magnetic fields of the neu¬ 
tron stars hosted in SFXTs (Grebenev & Sunyaev 2007; 
Bozzo et al. 2008), or a subsonic settling accretion regime 
combined with magnetic reconnections between the NS 
and the supergiant field transported by its wind (Shakura 
et al. 2012, 2014). 


^Project web page: http://www.ifc.inaf.it/sfxt/ 
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The Swift (Gehrels et al. 2004) 10 year anniversary 
also marks the completion of the first 7 years of the Swift 
SFXT Project, which has been investigating the properties 
of SFXTs with a strategy that exploits Swift's uniqueness. 
In this Paper I will review the main results of its long term 
monitoring and outburst follow-ups of an ever increasing 
sample of SFXTs, a monitoring that has recently extended 
to include a small control sample of classical systems. In 
particular, I will show the new results on the two sources 
IGR J16493-4348 and AX J1845.0-0433 which were ob¬ 
served during 2014, and put them in the broader context 
of the comparison of SFXT vs. classical systems. I will 
also show that, as customary, the seventh year harbors a 
“seventh year crisis”, since the most recent observations 
seem to challenge the prevailing models that attempt to 
explain the SFXT behaviour. 

2. Results from the Swift SFXT project 

A continuing effort to improve and fine-tune Swift's 
GRB observing strategy has allowed a gradual shift of 
overall observing time from mostly GRB science to (cur¬ 
rently) mainly guest observer and target of opportunity 
(ToO) targets. Several initiatives could then be carried 
out that would boost Swift's secondary science by select¬ 
ing well-defined astrophysical problems that could be effec¬ 
tively tackled by exploiting Swift's fast automatic slewing 
and multi-wavelength capability, as well as its flexible ob¬ 
serving schedule and very low overheads. The Swift SFXT 
Project was born as one of these initiatives. 
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Table 1: Summary of Swift/XJ{T campaigns divided by monitoring samples [yearly and orbital; Sect. 2.2 and 2.3). N is the number of 
observations (individual ObsIDs) obtained during the monitoring campaigns; Nc the number of observations used (Col. 6). Count rates 
(Col. 8) are in units of 10“^ counts s“^ in the 0.2-10 keV energy band, while observed fluxes (Col. 9) are in units of 10“^^ erg cm“^ s“^ and 
luminosities (Col. 10) in units of 10^"^ erg s“^, both in the 2-10keV energy band; all are based on a single 900s exposure. AT^; (Col. 11) is 
the sum of the exposures accumulated in all observations, each in excess of 900 s, where only a 3-cr upper limit was achieved; Pghort (Col. 12) 
is the percentage of time lost to short observations; IDC (Col. 13) is the inactivity duty cycle, the time each source spends undetected down 
to a flux limit of reported in Col. 9; RateAT^ (Col. 14) is the observed count rate in the data for which no detections were obtained as single 
observations (see Sect. 4). 


Name 

Period 

Distance 

Start 

End 

Nc/N 

Expo. 

lim 

lim 

r 2-10 

lim 

^short 

IDC“ 

ti^^..O'2-IO 

Rate^.^,^ 

Ref^ 

Ref'= 

Ref'^ 


(d) 

(kpc) 

UT 

UT 


(ks) 

(Xl0”3) 

(Xl0“^2) 

(X10^^) 

(ks) 

(%) 

(%) 

(X 10“^) 

P 

D 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 


Yearly sample 


IGR J08408-4503 

- 

3.4±0.35 

2011-10-20 

2012-08-05 

77/82 

74.4 

17 

1.9 

0.26 

46.6 

7 

67.2+J;? 

7.2ih0.6 

- 

12 

18 

IGR J16328-4726 

10.076 

6.5±3.5 

2011-10-20 

2013-10-24 

94/98 

88.0 

14 

2.7 

2.5 

47.5 

12 


4.0±0.4 

1 

13 

18 

IGR J16465-4507 

30.243 

12.7±1.3 

2013-01-20 

2013-09-01 

61/65 

58.6 

16 

2.0 

4.4 

3.0 

0 

5 1+4.4 
^•^-1.6 

14.6±0.4 

2 

12 

18 

IGR J16479-4514 

3.3193 

4.9 

2007-10-26 

2009-10-25 

139/144 

159.8 

16 

2.5 

1.1 

29.7 

3 

19 4+^'® 
^^•^-2.9 

3.1±0.5 

3 

14 

19,20,18 

IGR J16493-4348 

6.782 

> 6 

2014-01-19 

2014-09-03 

55/65 

52.7 

13 

2.2 

1.3 

5.8 

3 

11.3 tl'l 

> 7.9^ 

4 

15 

21 

XTE J1739-302 

51.47 

2.7 

2007-10-27 

2009-11-01 

181/184 

206.6 

13 

1.6 

0.18 

71.5 

10 

is.stll 

4.0±0.3 

5 

14 

19,20,18 

IGR J17544-2619 

4.926 

3.6 

2007-10-28 

2009-11-03 

138/142 

142.5 

12 

1.4 

0.21 

69.3 

10 

54.5+^;^ 

2.2ih0.2 

6 

14 

19,20,18 

AX J1841.0 — 0536 

6.4530 

7.8±0.74 

2007-10-26 

2008-11-15 

87/88 

96.5 

13 

1.8 

1.6 

26.6 

3 

28.4+|;i 

2.4±0.4 

7 

12 

19,18 
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5.7195 

6.4±0.76 

2014-02-14 
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13 

1.9 

0.11 

7.9 

4 

ii.at-« 

3.5±0.9 

8 

12 
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IGR J16418-4532 

3.73886 

13 

2011-02-18 

2011-07-30 

15/15 

43.3 

19 

12.5 

36 

4.8 

0 

> 9.2® 

9 

14 

22,18 

IGR J17354-3255 

8.448 

8.5 

2012-07-18 

2012-07-28 

22/22 

23.7 

14 

2.2 

3.3 

7.8 

1 

33 . 4 + 1 I 3 I 

> 4.6® 

10 

16 

23,18 

IGR J18483-0311 

18.545 

2.83±0.05 

2009-06-11 

2009-07-08 

23/23 

44.1 

11 

1.8 

0.24 

11.8 

0 

26.6 + 1“/ 

3.6 ± 0.8 

11 

17 

24,18 


^ Uncertainties obtained from or with the method described in Romano et al. (2014b). 

^ References to orbital periods: (1) Corbet et al. (2010); (2) La Parola et al. (2010); (3) Romano et al. (2009b); (4) Cusumano et al. (2010); (5) Drave et al. (2010); (6) 
Clark et al. (2009); (7) Gonzalez-Galan (2015) (8) Goossens et al. (2013); (9) Levine et al. (2011); (10) D’Al et al. (2011); (11) Levine &c Corbet (2006). 

^ References to distances: (12) Coleiro & Chaty (2013); (13) Piocchi et al. (2013); (14) Rahoui et al. (2008); (15) Nespoli et al. (2010); (16) Tomsick et al. (2009); (17) 
Torrejdn et al. (2010). 


d 


References to original data papers: (18) Romano et al. (2014a); (19) Romano et al. (2009b); (20) Romano et al. (2011a); (21) This work; (22) Romano et al. (2012); 
(23) Ducci et al. (2013b); (24) Romano et al. (2010). 


3-<7 upper limit. 






2.1. The long outburst of IGR J11215—5952 

The first experiment within the Swift SFXT Project 
was performed in February 2007 (Romano et al. 2007) 
when we took advantage of the fact that at least one 
SFXT, IGR J11215-5952 (hereon J11215), showed out¬ 
bursts with a periodicity of about 330 d (Sidoli et al. 2006) 
in the INTEGRAL data, and monitored one outburst with 
the S'm/t/X-ray Telescope (XRT, 0.2-10 keV, Burrows et al. 
2005) from before onset, to the peak, and until it be¬ 
came undetectable with reasonable XRT exposures. The 
campaign lasted 23 d (total exposure of ~ 73 ks) and 
showed that, differently from previously thought based on 
lower-sensitivity instruments observing only the brightest 
hr-long flares (e.g. the hard X-ray monitor on board IN¬ 
TEGRAL), the soft X-ray emission, hence the accretion 
phase, lasted several days. Superimposed on the light 
curve, frequent flares are also observed, probably due to in¬ 
homogeneous accretion, i.e. clumps in the accreting wind. 
This phenomenology is consistent with a gradually increas¬ 
ing flux at the periastron passage in a wide eccentric orbit. 

Swift also allowed the determination of the true orbital 
period Porb 164.6 d (Sidoli et al. 2007; Romano et al. 
2009c), a rare instance, as orbital periods are generally 
found from all-sky monitor data. 

Eurthermore, since the outbursts of this source are pe¬ 
riodic and regularly spaced, and a non negligible eccentric¬ 
ity is required to account for the observed low quiescent 
luminosity, these data led Sidoli et al. (2007) to propose 
the existence of a second wind component from the su¬ 
pergiant companion in addition to its normal symmetric 
polar wind. This equatorially enhanced wind component 
is probably clumpy, denser and slower than the polar one, 
and inclined with respect to the orbital plane (to account 
for the narrowness of the outburst X-ray light curve). 

2.2. Long term properties from monitoring campaigns 

In the wake of the success of the JII2I5 campaign, 
during the Fall 2007 we selected 4 SFXTs from the 8 
known at the time, IGR JI6479—4514, XTE J1739-302, 
IGR J17544-2619, and AX J1841.0-0536 (hereon J16479, 
J1739, J17544, and J1841, respectively, the seed of our 
yearly sample, see Table 1), and set to perform a system¬ 
atic study (Sidoli et al. 2008) with a year-long series of 
1-2 ks pointed observations with XRT. At the time, the 
SFXT binary periods were largely unknown and believed 
to be in excess of 10 d (see the current status in Table 1, 
Col. 2), so the observations were scheduled 3-4 days apart. 
The initial goals were to i) seek for the signatures of equa¬ 
torial winds, a) catch outbursts to determine whether they 
showed any periodicities. Hi) follow the sources during the 
whole outburst duration, and iv) monitor the quiescence. 

For this project, the Swift/R>nvst Alert Telescope (BAT, 
15-150 keV, Barthelmy et al. 2005) Team applied the “BAT 
special functions” to the sample of SFXTs and candidates 
(i.e., sources with similar X-ray flaring behaviour but with 
no firm measurement of the spectral type of the compan¬ 
ion). This allowed Swift to react to an increase in flux 


from any SFXTs as if they were a GRB. In this way, when¬ 
ever an SFXT triggers the BAT, simultaneous broad band 
data are collected that span ~ 1600-6000 A through the 
UV/Optical Telescope (UVOT, Roming et al. 2005), and 
0.2-150 keV through XRT and BAT combined. 

Before we began our investigation, deep XMM~Newton 
exposures (Gonzalez-Riestra et al. 2004) reporting fluxes 
between ~ 10“^^ and 10“^° erg had described the 
characteristics of J17544 away from the bright outbursts, 
including a trend for harder spectra at higher fluxes, while 
a revealing Chandra observation (in’t Zand 2005) had caught 
the first detection of a fast X-ray transient in quiescence, 
a state characterized by a very soft (photon index F = 
5.0 ± 1.2) spectrum. On the other hand, the long-term 
behaviour of SFXTs, not unlike any other hard X-ray 
transient, had been traditionally investigated only with 
coded-mask large field-of-view instruments, such as INTE¬ 
GRAL/IBIS (Ubertini et al. 2003) or Swift/BAT which, 
because of their sensitivity limits, mostly catch only the 
brightest portion of any transient event (see Romano et al. 
2014c, for a catalogue of more than a thousand Swift/BAT 
bright SFXT flares described in Sect. 2.5). Therefore, 
our strategy, by combining sensitive soft X-ray monitoring 
with outburst follow-ups, has allowed us through the years 
for the first time, to systematically assess the soft X-ray 
long term properties of a conspicuous fraction of the SFXT 
sample when away from the prominent bright outbursts, 
in particular, the states leading down to quiescence. 

The detailed results on the first campaigns can be found 
in Sidoli et al. (2008) and Romano et al. (2009b, 2011a), 
those on 3 more sources, IGR J08408—4503, IGR J16328 
-4726, and IGR J16465-4507 (hereon J08408, J16328, 
and J16465) in Romano et al. (2014a). Table 1 reports the 
characteristics of the sources in our SFXT sample (binary 
period and distance, Gols. 2 and 3) and the campaign dates 
(Cols. 4 and 5). In the following, we summarize the results 
for these sources and introduce the new ones for J16493 
and J1845, that were observed in 2014^. 

Fig. 1 shows the XRT light curves of the yearly sample 
at a daily resolution. The following can be observed: 

1. The most striking features are, of course, the bright 
outbursts, which will be described in detail in Sect. 
2.4. These outbursts, surprisingly, are however merely 
the tip of the iceberg in the SFXT activity (only a 
few percent of the total time). 

2. All sources display a dynamic range (DR) of ~ 3-4 
orders of magnitude^. 

3. The exceptions to the previous points are J16465 
and the newly observed J16493, the latter showing a 


^The new data were processed (see, e.g., Romano et al. 2014a) 
with standard software (FTOOLS v6.16), calibration (CALDB 
20140709), and methods (xrtpipeline, vO.13.1). 

^We note that both J1841 which went into outburst after the end 
of the campaign (2010 Jun 5, Romano et al. 2011b) and J1845 which 
went into outburst on 2012 May 05 (Romano et al. 2013), reached 
DR> 10^. 
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Figure 1: SwiftfKKT (0.2-10 keV) long term light curves for the yearly monitoring sample. The x-axis is in MJD. The points refer to the 
average count rate measured for each observation except for outbursts, where the data were binned to include at least 20 counts bin“^ (to 
best represent the dynamic range). Red points are detections, black downward-pointing arrows are 3-cr upper limits. Violet upward-pointing 
arrows mark bright outbursts detected by BAT (and not simultaneously followed by XRT), or by MAXI (for AX J1841.0—0536). The data 
on IGR J16493— 4348 (panel e) and AX J1845.0—0433 (panel i) are presented here for the first time. The remainder of the data are adapted 
from papers listed in Table 1 (Col. 17). 


DR~ 48. These two sources (Fig. Ih and le, respec¬ 
tively), based on their small DR, are therefore found 
to be not SFXTs, but classical systems. 

4. For SFXTs, most of the emission is found outside the 

bright outbursts, so that the long-term behaviour of 
SFXTs is not quiescence but an intermediate state 
with an average X-ray luminosity of erg s“^ 

(and, as reported below, a spectrum that can be fit 
with a power law with photon index F = 1-2). 

5. Variability is observed at all timescales we can probe. 
Superimposed on the day-to-day variability, we mea¬ 
sure intra-day flaring that involves flux variations up 
to one order of magnitude; we identify flares down 
to a count rate in the order of 0.1 counts s“^ {L ^ 2- 
6 X 10^^ erg s“^) within a snapshot of about Iks. 

As shown by Walter & Zurita Heras (2007) the short 
time scale variability cannot be accounted for by accre¬ 
tion from a homogeneous wind, but it can naturally be 


explained by the accretion of single clumps in the donor 
wind. We calculated that the average clump mass is Md ^ 
0.3-2 X 10^® g (Romano et al. 2011a), about those ex¬ 
pected (Walter & Zurita Heras 2007) to be responsible 
for short flares, below the INTEGRAL detection threshold 
and which, if frequent enough, may significantly contribute 
to the mass-loss rate. 

The data collected during these campaigns were used 
to perform soft X-ray intensity-selected spectroscopy. We 
find that the common out-of-outburst spectroscopic prop¬ 
erties are a non-thermal emission (power law with F = 
1-2) combined with a soft excess becoming increasingly 
more dominant as the source flux state becomes lower, and 
a ubiquitous harder-when-brighter trend (Romano et al. 
2009b, 2011a, 2014a). Therefore, the spectral modelling of 
out-of-outburst emission shows that accretion is occurring 
down to very low luminosities. This is clearly at odds with 
what is generally observed in the hard X-rays, due to the 
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Figure 2: Swift/WOT light curves of AX J1841.0—0536. Adapted 
from Romano et al. (2009b). 


different instrumental sensitivities involved. Since a non- 
thermal spectrum plus a soft excess is common in classical 
systems, the spectroscopic properties are not an efficient 
method for discriminating SFXTs within the HMXB sam¬ 
ple, differently from what happens for the dynamic range. 

We note that when combining single 1-2 ks snapshots 
in which no detection was individually achieved, the XRT 
data can reach luminosities comparable to the quiescent 
one (Romano et al. 2009b). For instance, in the cases 
of J17544 and J1739 (Romano et al. 2011a) and J08408 
(Romano et al. 2014a), luminosities of a few 10^^ ergs“^, 
have been reached, and the spectral properties observed in 
this very low state are consistent with those observed with 
deep XMM-Newton exposures (e.g. Bozzo et al. 2010). 

UVOT observed simultaneously with the XRT during 
most of our monitoring with different combinations of op¬ 
tical and UV filters, depending on the magnitude of the 
companion stars (Romano et al. 2009b, 2011a). In J1739, 
only marginal variability was observed in the u and uvwl 
filters; the uvwl light curve of J17544 and the u and uvwl 
light curves of J1841 were remarkably stable, consistently 
with the optical/UV emission being dominated by the con¬ 
stant contribution of the supergiant companions. Fig. 2 
shows, as an example, the UVOT light curves of J1841. 

2.3. Orbital monitoring campaigns 

Further monitoring campaigns were also performed on 
three more sources with higher-cadence pointed observa¬ 
tions (several XRT snapshots a day to provide intra-day 
sampling) for one or more orbital periods with the main 
goal of studying the effects of orbital parameters on the 
observed flare distributions. In 2009 we performed the 
very first complete monitoring of the soft X-ray activity 
along an entire orbital period (Porb ^ 18.5 d) of an SFXT, 
IGR J18483-0311 (J18483, Romano et al. 2010). Sim¬ 
ilarly, in 2011 we monitored the candidate SFXT IGR 
J16418—4532 (J16418, Romano et al. 2012) which has a 
much shorter orbital period, Porb 3.74d and, in 2012, 
the (then candidate) SFXT IGR J17354-3255 (J17354), 
with Porb = 8.4474d (Ducci et al. 2013b). These three 
sources compose our orbital monitoring sample. 

These unique datasets allowed us to constrain in these 
objects the different mechanisms proposed to explain their 
nature. In particular, we applied the clumpy wind model 



Energy (keV) 


Energy (keV) 


Figure 3: Swift/XIUt and BAT simultaneous spectroscopy. Filled 
(red) circles and empty (blue) circles denote XRT and BAT data, re¬ 
spectively. The data are fit with an absorbed power law with a high 
energy cut-off model or absorbed power laws with high energy expo¬ 
nential cut-off model. Data references: IGR. J08408—4503 (2008-07- 
05, Romano et al. 2009a); IGR J16328-4726 (2009-06-10, Romano 
et al. 2013); IGR J16479-4514 (2008-03-19, Romano et al. 2008); 
XTE J1739-302 (2008-08-13, Sidoli et al. 2009a); IGR J17544-2619 
(2009-06-06, Romano et al. 2009b); SAX J1818.6-1703 (2009-05-06, 
Sidoli et al. 2009b); AX J1841.0-0536 (2010-06-05, Romano et al. 
2011b); AX J1845.0-0433 (2012-05-05, Romano et al. 2013). 


for blue supergiants (Ducci et al. 2009) to the observed 
X-ray light curve. By assuming for J18483 an eccentricity 
of e = 0.4 and for J16418 circular orbits, we could ex¬ 
plain their X-ray emission in terms of the accretion from a 
spherically symmetric clumpy wind, composed of clumps 
with different masses, ranging from 10^® to xlO^^g for 
J18483, and from ^ 5 x 10^® g to 10^^ g for J16418. Since 
J18483 is an intermediate SFXT with a moderately high 
dynamic range in the X-ray luminosity, the estimated sizes 
and masses of the clumps in J18483 are somewhat larger 
that what would be expected according to the multidi¬ 
mensional simulations of massive stars winds (Dessart & 
Owocki 2002, 2003, 2005) but likely not unrealistic (see, 
e.g. Fiirst et al. 2014). The addition of magnetic/centrifu¬ 
gal gates could lower the requirements on the clump sizes 
and masses, but such mechanisms cannot be readily ap- 
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Figure 4: 5u;i/i/XRT light curves of the better followed-up out¬ 
bursts of confirmed SFXTs, referred to their respective BAT triggers 
(except: IGR J11215—5952 is referred to MJD 54139.94). Points 
mark detections, triangles 3a upper limits, vertical dashed lines mark 
daily Intervals, up to a week. Data references: IGR J08408—4503 
(2008-07-05, Romano et al. 2009a); IGR J11215-5952 (2007-02- 
09, Romano et al. 2007); IGR J16328-4726 (2009-06-10, Romano 
et al. 2013); IGR J16479-4514 (2005-08-30, Sldoli et al. 2008); 
XTE J1739-302 (2008-08-13, Sidoli et al. 2009a); IGR J17544-2619 
(2010-03-04, Romano et al. 2011b); SAX J1818.6-1703 (2009-05- 
06, Sidoli et al. 2009b); AX J1841.0-0536 (2010-06-05, Romano 
et al. 2011b); AX J1845.0-0433 (2012-05-05, Romano et al. 2013). 
Adapted from Romano et al. (2013). 


plied to the case of J18483 as the magnetic field of the 
NS hosted in this source and its spin period are highly 
debated (e.g. Diicci et al. 2013a; Sguera et al. 2015). We 
found that J17354 is probably a wind-fed system and that 
the dip observed in its light curve cannot be explained with 
a luminosity modulation due to a highly eccentric orbit; 
on the contrary, it can be explained in terms of an eclipse 
or the onset of gated mechanisms. 


2008). At the time of writing, we have collected a total 
of 51 bright flares (55 triggers, of which 4 double) that 
triggered the BAT, more than half of which followed-up 
with the NFl, thanks to the BAT special functions, so 
that we have been able to observe in this fashion most of 
the SFXT sample, as shown in Fig. 3. We found out that 
for all sources a good fit could indeed be obtained with a 
power-law with an exponential cutoff. This, in turn im¬ 
plies, based on the cutoff energy, a magnetic field consis¬ 
tent with a few 10^^ G, typical of accreting NSs in HMXBs. 

The availability of such data also motivated the devel¬ 
opment of a physical model, compmag in XSPEC (Farinelli 
et al. 2012a), which includes thermal and bulk Comp- 
tonization for cylindrical accretion onto a magnetized neu¬ 
tron star. A full description of the algorithm (see Farinelli 
et al. 2012a) is beyond the scope of this paper but the 
model has been successfully applied to the SFXT class pro¬ 
totypes J17391 and J17544 (Farinelli et al. 2012b), and to 
J18483 (Ducci et al. 2013a). 

An important benefit of NFl observations is Swiff's 
ability to provide arcsecond localization for several SFXTs 
and candidates whose coordinates were only known to the 
arcminute level, or to improve on previously known coor¬ 
dinates (e.g. Kennea et al. 2005; Grupe et al. 2009). This 
greatly helps in associating with optical counterparts. 

Our observing strategy also includes XRT follow-ups 
for days (generally up to a week) after the outburst via 
ToO observations, well after it had become undetectable 
with monitoring instruments with lower sensitivity. Fig. 4 
shows the best examples of outburst light curves as ob¬ 
served by XRT, and exemplifies the common X-ray char¬ 
acteristics of this class: 

• extended soft X-ray activity around an outburst last¬ 
ing up several days (see the vertical lines in Fig. 4); 

• a multiple-peaked structure; 

• a DR (only including bright outbursts) up to ^ 3 
orders of magnitude. 


2.4- SFXT broad-band properties and arcsecond localiza¬ 
tions from outbursts and outburst follow-ups 

Since Swift is a GRB-chasing mission, provided with 
fast automated slewing and panchromatic sensitivity, once 
SFXTs were included in the BAT special functions (see 
Sect. 2.2), SFXT outbursts started triggering the BAT 
and narrow field instrument (NFl) data (XRT and UVOT) 
started being collected within a few hundred seconds (down 
to about ~ 100 s) from the BAT trigger. The shape of the 
SFXT spectrum in outburst is a power law with an ex¬ 
ponential cutoff at a few keV, therefore the large Swift 
energy range can both help constrain the hard-X spectral 
properties (to compare with popular accreting neutron star 
models) and measure the absorption. 

These simultaneous XRT and BAT data allowed us to 
perform, for the first time, simultaneous broad band spec¬ 
troscopy of an SFXT in outburst (J16479, Romano et al. 


2.5. The 100-month SFXT BAT catalogue and the number 
of SFXTs in the Galaxy 

Since BAT observes an average of 88% of the sky daily, 
it is ideally suited to detect flaring hard X-ray astrophys- 
ical sources, SFXTs in particular. We have thus pro¬ 
duced the 100-month Swifi Gatalogue of SFXTs (Romano 
et al. 2014c) which collects over a thousand BAT flares 
from 11 SFXTs, and reaches down to 15-150keV fluxes 
of about 6 X 10“^° erg cm“^ s“^ (daily timescale) and 
about 1.5 X 10“® erg cm“^ s“^ {Swift orbital timescale). 
We found that these hard X-ray flares typically last at 
least a few hundred seconds, reach above lOOmCrab (15- 
50 keV), and last much less than a day. Their clustering 
in the binary orbital phase-space, however, demonstrates 
that these short flares are part of much longer outbursts, 
lasting up to a few days, as previously observed during 
our outburst follow-ups (Sect. 2.4). This large dataset can 
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Figure 5: Differential distributions of the 2-10 keV luminosity (lower axis) and flux (unabsorbed, upper axis) drawn from the Swift/X.Kr 
light curves binned at 100 s. N is the sample size. The thick black lines represent the data collected during the monitoring campaigns (see 
Table 1, Col 4 and 5 for the campaign dates). The thin blue histograms show the part of the data collected as outburst observations during 
the campaign, thus including both the initial bright flare and the follow-up observations (see Sect. 3). The thin red histograms show outburst 
observations collected outside of the monitoring campaign (one outburst per source). The data on IGR J16493—4348 and AX J1845.0—0433 
are presented here for the first time, the rest are adapted from Romano et al. (2014a). 


therefore probe the high and intermediate emission states 
in SFXTs, and help infer the properties of these binaries; it 
can also be used to estimate the number of flares per year 
each source is likely to produce as a function of the detec¬ 
tion threshold and limiting flux in futnre missions. Finally, 
the catalogue has recently been exploited by Ducci et al. 
(2014) to estimate the expected number of SFXTs in the 
Milky Way, N ^ 37^22- This shows that SFXTs con- 
stitnte a large portion of X-ray binaries with supergiant 
companions in the Galaxy. 

3. Differential luminosity distributions 

An often understated property of the monitoring data 
is that the yearly campaigns are statistically representa¬ 
tive of the long-term soft X-ray properties of SFXTs that 
the deep exposures from pointed telescopes can only rarely 
and non-uniformly sample. Our observations are also in- 
dependenf since each observation is not triggered by the 
preceding ones (we consider the outburst followups as a 
separate set of data when in need of a statistical set). Onr 


monitoring pace thus provides a casual sampling of the soft 
X-ray light curve at a resolution of ^ 3-4 d over a base¬ 
line of one or two years, therefore it ofl’ers both coverage 
of a large number of binary orbital cycles (ranging from 
^ 15 cycles for J17391 to ~ 220 for J164794) and a good 
sampling of the orbital phase (see fig. 6 of Romano et al. 
2014a). 

Based on these premises, we can effectively calculate 
the percentage of time each source spent in different flux 
states, among which we distingnish: 

1. the flares that trigger the BAT (see Sect. 2.4) ac¬ 
counting for 3-5% of the exposure time; 

2. the intermediate states (all observations yielding a 
firm detection exclnding ontbursts); 

3. non detections (significance below 3a; see Sect. 4). 

Only observations with an exposure in excess of 900 s 
were considered to account for non detections ob¬ 
tained during very short exposnres (dne to our obser¬ 
vations being interrnpted by a higher figure-of-merit 
GRB). These correspond to flux limits ^ 

(1-3) X 10“^^ erg cm“^ s“^ (see Gol. 9 of Table 1, and 
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also the corresponding limits in count rates, Col. 8, 
and luminosity, Col. 10). 

From the XRT light curves binned at 100 s, of both 
the yearly and the orbital monitoring samples, and af¬ 
ter removing the observations where a detection was not 
achieved (thus selecting only the intermediate states), we 
construct the 2-10 keV differential luminosity distributions 
(DLD), shown in Fig. 5 (solid black lines). For all sources 
we adopted a single conversion factor between count rates, 
fluxes and luminosities, that were derived from the ‘medium’ 
spectrum for J16465, J16479, J16493, J1739, J17544, J1841, 
J1845, and J18483, the ‘low’ spectrum for J08408 and 
J16328, the first observation for J16418, and the average 
spectrum for the weak source J17354 (Romano et al. 2010, 
2011a, 2012; Ducci et al. 2013b; Romano et al. 2014a). Be¬ 
cause the uncertainty in this conversion is dominated by 
those on the distance (Table 1, Col. 3), the top x-axis of 
Fig. 5 also reports the flux scale in the same energy band. 
We note that the DLDs drawn from the orbital monitoring 
sample need to be taken with caution, as they follow an 
entirely different observing strategy. These observations 
were in fact collected with intensive campaigns during one 
or a few orbital periods, as opposed to the few points per 
each period typical of the yearly campaign data, so short 
timescale variability may play an important role. 

Fig. 5 distinguishes (as a thin blue line) the data that 
were taken during an outburst that occurred during the 
observing campaign (two for J16479, and three for J1739 
and J17544), and those (thin red line) of outbursts that 
were observed outside of the campaigns (we considered one 
for each of J08408, J16328, J1840, and J1845). 

The DLD of the SFXT prototypes, J1739 and J17544, 
as well as J16479 and J08408 show two distinct popula¬ 
tions of flares. The first one is due to the outburst emission 
(peaking/reaching a few 10“® erg cm“^ s“^), the second 
is due to the out-of-outburst emission, characterized by 
emission spanning up to 4 orders of magnitude in DR (at 
100s binning). This also applies to the newly observed 
J18450, which also shares a DR of at least 3 orders of 
magnitude at a temporal resolution of 100 s. We cannot ex¬ 
clude that particular distributions of the clump and wind 
parameters may produce a double-peaked DLD (Romano 
et al. 2014a), but this behaviour is more easily explained 
in terms of different accretion regimes as predicted by the 
magnetic/centrifugal gating model or the quasi-spherical 
settling accretion model (Grebenev & Sunyaev 2007; Bozzo 
et al. 2008; Shakura et al. 2012, 2013). 

The classical systems (J16465 and the newly observed 
J16493), on the contrary, only show one peak in their 
DCD, which is significantly brighter than those of of SFXTs. 
This confirms the findings of Lutovinov et al. (2013) that 
SFXTs show a median luminosity beneath the one of nor¬ 
mal wind-fed HMXBs; the flaring observed in SFXTs can 
therefore be explained if some mechanism, such as mag¬ 
netic arrest, can inhibit accretion. DCDs, therefore, can be 
used effectively to discriminate between the most extreme 
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Figure 6: Distribution of the XRT detections (count rates) of 
IGR J17544—2619 (0.2-10 keV) folded at the orbital period, based 
on the most recent outburst ephemeris for this source {P = 4.92693di 
0.00036 d, periastron at MJD 53732.65 zb 0.23; Smith 2014). 

SFXTs, the intermediate systems, and classical systems. 

While the outbursts account for 3 and 5 % of the expo¬ 
sure time, the most probable flux level at which a random 
observation will find these sources, when detected, can be 
retrieved from the peak of the DLD, and this is of course 
a powerful tool to plan further observing campaigns, as is 
the distribution of detections along the orbital phase (see 
fig. 6 of Romano et al. 2014a). In particular, in Fig. 6 
we show the case of J17544 which flares tend to cluster 
around periastron more than in other SFXTs. 

4. Inactivity duty cycles 

The duty cycle of astrophysical sources is usually de¬ 
fined as the fraction of time the sources are active, and it 
is used to both characterize their emission properties and 
and to plan further observing campaigns to study them. 
SFXTs, however, show a very large dynamical range, with 
activity observed by the XRT spanning several orders of 
magnitude in flux. It is more interesting, therefore, to 
define a measurement of inactivity as opposed of one of 
activity. 

From the non detections, we dehne the inactivity duty 
cycle (Romano et al. 2009b) as the time each source spends 
undetected down to a flux limit of 1-3 x 10“^^ erg cm“^ s“^, 

IDC = Ars/[ATtot (1 - Pshort)] , (1) 

where ATs is the sum of the exposures (each longer than 
900 s) accumulated in all observations where only a 3 (t 
upper limit was achieved (Table 1, Col. 11), ATtot is the 
total exposure accumulated (Table 1, Col. 7), and Pshort 
is the fraction of time lost to short observations (exposure 
< 900s, Table 1, Col. 12). The cumulative count rate for 
each object is also reported Table 1 (Col. 14). 

The need to provide uncertainties on IDCs by avoiding 
the standard approach of deriving them from extensive 
and time-consuming Monte Carlo bootstrap simulations 
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Figure 7: The XRT duty cycle (2-10 keV) as a function of orbital 
period and for a range of 2-10 keV luminosities (see legend, in units 
of erg s“^) in black. Only points above 0.1% are shown. Since the 
period of J08408 is currently unknown (see Table 1 Col. 2 for the 
values adopted for the rest of the sample), we arbitrarily place it at 
an orbital period of about two days and mark it in blue. The shaded 
areas mark the loci of the XRT duty cycle defined by contiguous 
luminosities. The red filled stars represent the XRT duty cycle at 
the INTEGRAL sensitivity for each object (the downward-pointing 
arrows are consistent with 0). The data on IGR J16493—4348 and 
AX J1845.0—0433 are presented here for the first time, the rest are 
adapted from Romano et al. (2014a). 

has lead us to propose an application of Bayesian tech¬ 
niques, instead (Romano et al. 2014b). We exploited the 
fact that SFXTs are, when considering duty cycles, two- 
state sources, since they can only be found in either of 
two possible, mutually exclusive states, inactive (off) or 
flaring (on). Or, in this case, above or below a given flux 
threshold. We derived the theoretical expectation value for 
the duty cycle and its error as based on a finite set of in¬ 
dependent observational data points following a Bayesian 
approach (Romano et al. 2014b). The IDCs and their un¬ 
certainties thus calculated are reported in Table 1, Col. 13. 

We note that IDCs can be quite small for classical 
systems, which are to all extents and purposes persistent 
sources, and because they are on average more luminous. 
Once again, IDCs can help discriminate between SFXTs, 
intermediate and classical systems. 

5. Discussion: the seventh year crisis 

In the following, we discuss the “seventh year crisis”, 
the challenges that the recent observations (those collected 
by Swift, in primis) are making to the prevailing models 
attempting to explain the SFXT behaviour. 

Duty Cycles and orbital geometry. If the proper¬ 
ties of the binary geometry and inhomogeneities of the stel¬ 
lar wind from the primary were the leading causes of the 
observed X-ray variability in SFXTs, as initially proposed 
in clumpy wind models (e.g. in’t Zand 2005; Negueruela 
et al. 2008; Walter & Zurita Heras 2007), generally larger 


IDCs would be expected for longer orbital periods. Natu¬ 
rally, the definition of duty cycle is strongly dependent on 
the luminosity assumed as lower limit for the calculation. 
Therefore, we exploited the high sensitivity afforded by the 
XRT observations and defined an XRT luminosity-based 
duty cycle (XRTDC) as the percentage of time the source 
spends above a given luminosity. We considered several 
(2-10 keV) luminosities in the range L 2 _iokeV = lO^^-IO^® 
erg s“^ and included the particular value of the luminos¬ 
ity corresponding to the INTEGRAL sensitivity for each 
object. 

Figure 7 shows the XRTDC as a function of the orbital 
period with red stars marking the value at the INTEGRAL 
sensitivity. No clear correlation is found between the or¬ 
bital periods and any of the duty cycles. This implies 
that wide orbits are not characterized by low duty cycles, 
as the clumpy wind models would predict. An intrinsic 
mechanism instead seems to be more likely responsible for 
the observed variability in SFXTs, i.e., either the wind 
properties or the compact object properties. However, it 
is hard to justify the radically different wind properties 
in SFXTs from those in classical systems with the same 
companion spectral type. Therefore, in light of this lack 
of correlation with the orbital period and our finding dis¬ 
tinct flare populations (see Sect. 3), it seems more plausible 
that accretion-inhibition mechanisms or a quasi-spherical 
settling accretion regime may be in action, instead. 

Cumulative luminosity distributions. In Sect. 3 
and 4 we have shown that DLDs and IDGs can be used 
effectively to discriminate SFXTs from classical systems, 
since the former are characterized by lower average lu¬ 
minosities. Another way to examine this property is to 
use the cumulative luminosity distributions (GLD) for our 
sample as calculated from the long term monitoring data 
(Sect. 3). In Bozzo et al. (2015) the GLDs of the SFXT 
sample, as well as that of the classical SGXB IGR JI8027 
—2016 (JI8027), are reported in the soft X-rays. Previ¬ 
ous work constructed GLDs based on the RXTE Galactic 
bulge scan programme data (Smith et al. 2012) and IN¬ 
TEGRAL long-term monitoring (Paizis & Sidoli 2014). In 
Fig. 8, in which we also added the newly observed JI6493 
and JI845, the GLDs are normalized to the total exposure 
for each source, so that the source duty cycle corresponds 
to the highest value on the y-axis. We can see that classical 
systems are characterized by GLDs with a single knee at 
^ lO^^-lO^"^ erg s“^. On the contrary, SFXTs are system¬ 
atically sub-luminous, with their GLDs shifted at lOO-IOO 
times lower luminosities. 

As shown in Fig. 8, the classical SgXRB JI8027 (thick 
dashed line) is characterized by a single knee at ~ 10^® 
erg s“^. The GLDs of JI6465, JI64I8, and JI6493 (thick 
dashed lines) closely resemble that of classical systems, 
with a position of the knee that can be accounted for once 
the relative distances and dependence of the flux from the 
orbital periods are considered. As done in Bozzo et al. 
(2015) for the former two, we here reclassify the newly 
observed SFXT candidate JI6493 as a classical system. 
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Figure 8: Cumulative luminosity distributions of the 2—lOkeV luminosity drawn from the Swift/'X.Kr light curves binned at 100 s. The 
classical SGXB IGR 18027—2016 is marked with a thick dashed line. Newly classified classical systems IGR J16465—4507, IGR J16418—4532, 
and IGR J16493—4348 are also marked with thick dashed lines. The intermediate SFXTs IGR J18483—0311 and IGR J17354—3255 are shown 
as dot-dashed lines. IGR J16328—4726, IGR J16479—4514, and AX J1841.0—0536 are shown as dotted lines, while the most extreme SFXTs 
(IGR J08408—4503, XTE J1739—302, IGR J17544—2619, and AX J1845.0—0433) are marked with solid lines. The data on IGR J16493—4348 
and AX J1845.0—0433 are presented here for the first time, the rest are adapted from Bozzo et al. (2015). Left: Only XRT data collected 
during the monitoring campaigns of all sources are used. Right: Same as for the left, but in this case we also considered one outburst for 
the sources IGR J08408—4503, IGR J16328—4726, AX J1841.0—0536, and AX J1845.0—0433 recorded by the XRT outside the corresponding 
monitoring campaigns. 


The intermediate SFXTs J18483 and J17354 (dot-dashed 
lines) have CLDs similar to those of classical systems, 
but shifted at ten times lower luminosities. The CLDs of 
J16328, J16479, and J1841.0 (dotted lines), are character¬ 
ized by even lower luminosities (one can note the similarity 
of the orbital periods of J16479 and J16418, and that of 
J16328 to Vela X-1). The profiles show more complexity, 
as more knees are appearing, reflecting different peaks in 
the DLDs, that is, different population of flares. Finally, 
the most extreme SFXTs, J08408, J1739, J17544, and the 
newly observed J1845.0 (solid lines) show very complex 
profiles. 

Considering that both differential/cumulative distribu¬ 
tions and inactivity duty cycles show that SFXTs are un- 
derluminous when compared to HMXBs, and that single 
knee profiles in CLDs can be understood in terms of wind 
accretion from an inhomogeneous medium (see, e.g. Fiirst 
et al. 2010), we can interpret the differences observed be¬ 
tween classical systems and SFXTs as due to accretion 
from a structured wind in the former sources and to the 
presence of magnetic/centrifugal gates or a quasi-spherical 
settling accretion regime in the latter. 

The king, the power and the ring. Ruler of 
a small kingdom, the SFXT prototype J17544 has been 
a stimulus, a catalyst, and a continuous challenge in the 
process of understanding SFXTs as a class since it was 
discovered by INTEGRAL over ten years ago (Sunyaev 
et al. 2003). The optical counterpart in this binary (or¬ 


bital period Porb = 4.926 ± 0.001 d, Clark et al. 2009) is 
quite an ordinary 09Ib star with a mass of 25-28 Mq (Pel- 
lizza et al. 2006) and located at a distance of 3.6 kpc (Ra- 
houi et al. 2008). The large luminosity swings observed on 
timescales as short as hours were alternatively explained 
by mechanisms that regulate or inhibit accretion (Stella 
et al. 1986; Grebenev & Sunyaev 2007, propeller effect; 
Bozzo et al. 2008, magnetic gating). In particular, Bozzo 
et al. (2008) explained them in terms of transitions across 
the magnetic and/or centrifugal barriers. In this scenario, 
the large dynamic range of SFXTs (about five decades) 
is achieved with a small variation of the mass loss rate 
(e.g. a factor of ^5 in fig. 3a of Bozzo et al. 2008) by as¬ 
suming a spin period Pspin in excess of ^ 1000 s and a 
magnetar-like field (P > lO^'^G). A recent NuSTAR ob¬ 
servation (Bhalerao et al. 2015), however, has revealed a 
cyclotron line at 17keV, yielding the first measurement of 
the magnetic field in a SFXT, at ^ 1.5 x 10^^ G, typical of 
accreting NSs in HMXBs. This set of observations there¬ 
fore rule out the magnetar nature for the SFXT prototype. 

On the other hand, the propeller and gating models are 
still applicable to fast rotators (Pspin < 10^ s) even with 
non-magnetar magnetic fields like those observed in classi¬ 
cal systems. Grebenev (2010) provided a simple equation 
to estimate the neutron star spin period at which the mag¬ 
netic inhibition regime takes place. 
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Figure 9: Bolometric X-ray luminosity light curves of the bright¬ 
est outbursts recorded by Swift for IGR J17544—2619. The giant 
burst of 2014 October 10 is shown in red, and compared with pre¬ 
vious bright outbursts from this source. The horizontal (pink) line 
marks the average level for this source, obtained from the two year 
monitoring campaign (Romano et al. 2011a, see Sect. 2.2), while the 
right-hand y-axis is the standard Eddington luminosity for spherical 
accretion of fully ionized hydrogen for a 1.4 M0 NS. Adapted from 
Romano et al. (2015). 



where B 12 = B/10^^ G is the magnetic field of the neutron 
star, M _5 = Mq yr“^ and V 3 = Tw/ 10^ km 

s“^ are the mass loss rate and the wind velocity of the 
donor star, and Pio = Poih/W d is the orbital period (see 
also fig. 4 of Grebenev 2010). Bozzo et al. (2008) showed 
that for sub-magnetar fields and < 100 s spinning neutron 
stars, a 10® luminosity swing would require a much higher 
increase in the mass loss rate with respect to the magnetar 
case (a factor of ^ 100 in fig. 4 of Bozzo et al. 2008). 
Such increase would imply that, either the mass loss rate 
from the supergiant stars in SFXTs changes abruptly on 
short time scales, or the local fluctuations in velocity and 
density of their winds are substantially larger than that of 
classical systems. In both cases, this would imply a fine 
tuning for the properties of the supergiant star winds in 
SFXTs compared to other SGXBs. At present, there is 
no observational evidence that this should be the case, as 
discussed in Bozzo et al. (2015). 

While current investigations concentrate on finding pos¬ 
sible mechanisms to inhibit accretion in SFXTs and to ex¬ 
plain their unusually low average X-ray luminosity, J17544 
seems to be testing our modelling further. An exception¬ 
ally bright outburst was observed by Swift on 2014 October 
10 (Romano et al. 2015), during which the source reached 
a peak luminosity of 3 x 10®®ergs“^ (or a 0.3-10keV un¬ 
absorbed flux of 1.0 X 10“^ erg cm“^ s“^, corresponding 
to 2.1 Grab). Tentative evidence for pulsations at a period 
of 11.6 s, and an expanding X-ray halo, or ring, around 
the source were also found in the XRT data. Such a high 
luminosity (see Fig. 9) not only extends the dynamic range 
of this source to DR^ 10®, a uniquely high value (by a fac¬ 
tor of 10), but also reaches the standard Eddington limit 
expected for a NS of 1.4 Mq, thus challenging, for the first 


time, the maximum (rather than the minimum) theoret¬ 
ical luminosity expected for an SFXT. In Romano et al. 
(2015) we propose that this giant outburst could be caused 
by the formation of a transient accretion disc around the 
compact object. 

Conclusions 

In the last seven years. Swift has contributed many 
“firsts” in the SFXT field, that I have summarized in this 
Paper; most importantly, it has performed the first sys¬ 
tematic investigation of the soft X-ray long term prop¬ 
erties of SFXTs with a vey sensitive instrument. This 
has provided us with many clues to help us understand 
SFXT outburst physics; it has revised or revolutionized 
incomplete or over-inferred properties derived from lower- 
sensitivity monitorings; it has motivated the creation of 
new models (both geometrical and physical) and helped 
test their applicability. 

Swift has consistently surprised us with the unexpected. 
And yet, we are still missing some key ingredients to un¬ 
derstand SFXT variability, especially when compared with 
classical systems. The current crisis is an excellent moti¬ 
vation to look “deeper and longer” and, in this frame¬ 
work, Swift monitoring programs on SFXTs and classical 
HMXBs will be crucial. 
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